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Ch ■ Abstract 

We show that the cosmological phase transition from the first accelerated expansion in the 
I early universe to the second accelerated expansion over the intermediate decelerated expansion 

is possible in the HL gravity without the "detailed balance" condition if the dark scalar energy 
pi • density is assumed to be negative. Moreover, we obtain various evolutions depending on the scale 

, factor and the expansion rate. Finally, we discuss the existence of the minimum scale in connection 

_ with the singularity free condition. 
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I. INTRODUCTION 



Recently, the Horava-Lifshitz (HL) gravity has been proposed as an ultraviolet (UV) 
completion of general relativity motivated by the Lifshitz theory in the condensed matter 
physics The key of the UV completion is the anisotropic scaling between space and time, 

t^¥t, x'^bx\ (1) 

where the Lifshitz parameter z becomes 1 in the infrared (IR) limit. In general, the HL 
theory is not invariant under the full diffeomorphism group of general relativity but under 
its subgroup, called the foliation-preserving diffeomorphism; however, in the IR limit, the 
full diffeomorphism is somehow recovered. A mechanism for recovering the full diffeomor- 
phism or the renormalization group flow is yet unsolved issue, but the HL gravity has been 



intensively studied in the area of cosmology |3l-lll| and black hole physics 12l-ll9|. 



Considering Arnowitt-Deser-Misner (ADM) decomposition of the metric with ds 
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-N'^c^dt^ + gij{dx' + N'dt){dx^ + N^dt) |20|, the Einstein-Hilbert action can be rewritten 
as ^ 

Ieh = T77^ / dS^ [n - 2A] 
c2 r 

= TTTV^ / dtd^^^N [K,jK'^ - + c\R ~ 2A)] , 
IbTrGAT J 

where Kij = [gij — ViNj — VjNi] is the extrinsic curvature of t = constant hyper- surface, 
and the dot denotes the derivative with respect to t. Here, Qij, R, and Vj are the metric, 
the intrinsic curvature, and the covariant derivative in the three-dimensional hyper-surface, 
respectively. Note that the action (jl]) can be regarded as one of the z = 1 HL theory and 
each terms are invariant under the foliation-preserving diffeomorphism. Then, the first two 
terms are referred to as kinetic terms, while the other two are potential terms. 

In particular, to study the HL gravity, higher-order potential terms, for instance, R^, 
will be taken into account. But there are almost ten possible terms for z = 3 case, and 
the "detailed balance" condition can reduce the ten coefficients to three effective ones 
Of course, it also implies that {D + l)-dimensional renormalization can be reduced to the 
simpler dimensional renormalization. By the way, it has been claimed that matter is not 
UV stable with this condition [2l|, and the HL theory with the detailed balance condition 
does not have the Minkowski vacuum solution. 
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On the other hand, there have been interesting studies on the cosmological evolution in 
the HL gravity based on the phase space analysis jl^. Authors showed that the cosmolog- 
ical phase changing from the decelerated expansion to the accelerated expansion is possi- 
ble assuming the detailed balance condition and the constant equation-of-state parameter. 
However, the first accelerated expansion corresponding to the inflationary era has not been 
discussed. We believe that there should be such a structure if the HL theory is indeed the 
UV completion. By the wa y, t he HL gravity itself may be problematic, for instance, counting 
its degrees of freedom 22, 23], description of the asymptotically flat spacetime p, Q, and 



so on. However, even in spite of these problems, it deserves to study various aspects since it 
may give some insight to understand the quantum gravity. So, we would like to investigate 
whether the smooth cosmological phase transition from the first accelerated expansion in 
the early universe to the second accelerated expansion over the intermediate decelerated 
expansion is possible or not. 

In section [Tll we recapitulate the Hofava-Lifshitz cosmology and define some relevant 
quantities without the "detailed balance" condition. Then, the nonlinear equations of mo- 
tion will be solved using nonlinear methods in order to study the cosmological behavior 
in section lllli By considering the observational data on density parameters, we shall plot 
the phase portrait in section IIVI and discuss the behavior in the early stage of our uni- 
verse and its various destiny. If a dark scalar energy density is assumed to be negative, 
then the desired first accelerated expansion appears. Moreover, we will show that it can 
be smoothly connected with the second accelerated expansion. Finally, discussions will be 
given in section |Vl 



II. COSMOLOGICAL SETTING FOR z = 3 

We now consider the anisotropic scaling ([T]) between time and space, and then gij and 
are invariant while A^* — ?■ b^'^N^ and c — t- b^^^c. Requiring that the Planck constant is 
also invariant, h ^ h, we get E — )■ b^^E and M — )■ b^~'^M, where E and M are energy and 
mass, respectively. Then, the kinetic action in the Horava-Lifshitz gravity (HL) is naturally 
given as 

/kin = 4 / dtSx^N [K,^K'^ - XK^] , (3) 
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where is a coupling related to the Newton constant Gn-, and A is an additional dimen- 
sionless coupling constant. Note that the original kinetic part of Einstein- Hilbert action ([2]) 
can be recovered when A = 1 and = 32'7rGAr/c^. Moreover, the power-counting renor- 
malizability requires z > 3 [l^. From now on, we will choose 2; = 3 for simplicity, then 
becomes dimensionless because of /t^ 



Now, the most general potential for z = 3 can be written as 



'pot 



— / dtSx^Nya + (5R + 71 + ^^R^^R^ + ^e''''RaV,Ri 

+ aiR^ + a2RRiiR'^ + a^RlR^R^ + a^S/iRV'R + a^V iR^jk^' R^^ 



(4) 



Note that other possible terms are not independent because of the Bianchi identity and 
symmetries of Riemann tensors. The potential terms in Eq. ([2]) can be recovered when we 



set a = 2Ac^ (3 
identified with 



-c^ and 'Ji = ^ = (Ti 



— /3, Gn 



0. Then, the fundamental constants can be 



32^'' 



A 



a 



(5) 



On the other hand, the ten coefficients can be reduced to three, /x, A^, (, using the 
detailed balance condition 



a 



3fi:V^AH/ _ /tV^AvF _ - 4A) 



16(3A - 1) 



16(3A - 1^ 



7i 



^4 



cr2 



5^4 



64(3A - 1) 
3k^ 



72 



4 

K H 
■4^' 



K 



(6) 



■32C 4(4' 



for which the three-dimensional topologically massive gravity action [24,] is given by 



W = fi J d'x^{R-2Aw) + ^ J 



x(r), 



(7) 



where x(r) represents the gravitational Chern-Simons term. However, we will take the 
general ten coefficients without resort to the detailed balance condition (|6]) so that the total 
action becomes 

I = Ikin + -^pot + Imat (8) 

where we introduced a normal matter action /mat which is a perfect fluid source of the energy 
density p and the pressure p. 

Now, we are going to consider the Robertson- Walker (RW) metric. 



ds' 



-c^dr + a\t) 



+ r^dni 



(9) 
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where k = 0, ±1 is the normahzed spatial curvature. After some tedious calculations in 
Eqs. dH]), the equations of motion can be obtained as 

2 _ ^2 , ^ [a , k(3 , 2A;2(37i + 72) , 4:k'\9ai + 3(T2 + as] 



3{3\-l)H' = —p + 6 



6 ^ a2 ^ a4 ^ a« 



Y [P + Pvac + Pk+ Pdv + Pds] , (10) 
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A;/3 , 2^2(871 + 72) , 4A;3(9ai + 3(72 + (73) 



+ — — + 



6 3a2 3a4 



(3A-1) (^+^^^)= 

= - Y + ^^'^c + + ^'i'' + , (11) 

where H = d/a is the Hubble parameter. Apart from the normal matter contribution, 
the additional energy-momentum contributions from the potential terms (jl]) are explicitly 
written as 

_ 3c2 a _ 1 _ c2 

Pvac Pvac 7:; "T) Pk T^Pk 7:; 

_1 c2 2A;2(37i+72) _ 3c2 4^=^(9^1 + 3^2 + ^3) 

Pdr - gPdr - g^^^ ^4 , Pds - Pds - g^^^ ^6 

Similarly to the general relativity, p^ac and pk come from the vacuum energy and the spatial 
curvature contributions, respectively. In particular, in the HL theory, pdr is called the dark 
radiation. Moreover, pds is a dark scalar, since it is characterized by pds = Pds ~ It 
is interesting to note that for the spatial flat geometry, the equations of motion (fTOj) and 
f fTT]) are reduced to those of general relativity up to a factor (3A — l)/2; in other words, the 
given HL cosmology is prominent for the nonvanishing spatial curvature. Both pdr and pds 
depend on the spatial curvature k; in this sense, these might be regarded as corrections to the 
curvature contribution pk- The total actual energy can be defined by ptot = P + Pvac + Pdr + Pds- 
From Eqs. ( |T0|) and ( ITTl) . we can get the acceleration of the scale factor, 

(3A-l)^ = -^[ptot + 3ptot], (13) 

where we used the relation p^ + 3pfc = 0. Next, differentiating Eq. f IlOp with respect to t and 
plugging it into Eq. (fTTj) . we can naturally obtain the fluid equation pc + 3H {pc + Pc) = 0, 
where pc is defined by pc = 6(3A — 1)H'^ / = ptot + Pk- Without loss of generality, the 
conservation equation becomes 

Pi + 3H {pi + Pi) = 0, (14) 
which is valid for each source component, i.e., {pi} = {p, Pvac, Pfc, Pdr, Pds}- 



III. LINEAR ANALYSES 



Ultimately, we want to show the cosmological phase transition from the first accelerated 
expansion to the second accelerated expansion in terms of the intermediate decelerated 
expansion, which looks like the overall behavior of our universe. Before we get down to this 
problem, we are going to exhibit the essence of the linear analysis. First, rewriting Eq. (fT3|l . 
we have 

2 

(3A-l)^ = -^[l + 3a;(a)]ptot, (15) 

where a; (a) = ptot/ptot is the equation-of-state parameter for the total energy. To describe 
phase portraits, we should make change of variables so that we introduce v = a, and h = 
K^ptot/Q{SX - 1) with A > 1/3. Then, we can get three first-order differential equations. 



a = f , 



V = — - [1 + 3ijj{a)] ah, 



(16) 



h = -3 [l + uj{a)] 



vh 



where we used Eqs. f lT^ and f|T5l) . Note that integrating Eq. ([T6l) . we can find a conservation 
relation, 

2 I, 2 2k/3 f-irrS 

ah — v= —— r, (17) 



3A- 1' 

where the integration constant in the right hand side of ( |T71) was consistently fixed by 
comparing Eq. (|T0|) . Now, let us eliminate h in Eq. f|T6|) using the conservation relation ([T7|) . 
then we have just two first-order differential equations. 



a = v = f{a,v), 



g[a,v) 



(18) 
(19) 



where the equation-of-state parameter can be written as u{a) = — 1 — ah' /3h with h' = 
dh/da. 

To study the behavior of the scale factor in the phase space, we should calculate the 
Jacobian matrix. 



A = 


daf d^f 








dag d^g 


{a*,v*) 



1 

3kl3uj'{a*)/{3\ - l)a* 



(20) 
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zero. As for k ^ 0, the Jacobian matrix (12 Op becomes A 



, and 



on the fixed point {a*,v*) which is easily obtained from a = -y = in Eqs. (fTSjl and (fT9|) . 

f V* = 0, for k = 0, 

(a*, V*) satisfying < (21) 

[l + 3uj{a*) = 0,v* = 0, hrk^O. 

Note that for k = 0, there are non-isolated fixed points; i.e. every point on the line v* = 
can be a fixed point. Next, considering an eigenvalue equation of = 77V which gives 
two eigenvalues, the fixed points can be classified as follows 25|]: They are repellers if both 
eigenvalues have positive real part, attractors if both eigenvalues have negative real part. In 
particular, they become saddles if one eigenvalue is positive and the other is negative. On 
the other hand, if both eigenvalues are pure imaginary, then they are centers. Moreover, 
they are called non-isolated fixed points if at least one eigenvalue is zero. 

Specifically, the Jacobian matrix fl20|) for A; = 0, is trivial so that we can find non- 
isolated fixed points as mentioned earlier, due to the fact that at least one eigenvalue is 

1 

3A;/3w'(a*)/(3A - l)a* 
the corresponding eigenvalues are 77 = )/(3A — l)a* where A > 1/3 is already 

assumed. Then, the fixed point ( 12T1) becomes a saddle node for k(3u'{a*) > or a center for 

k(3uj'{a*) < 0. 

IV. COSMOLOGICAL PHASE TRANSITION 

We are now in a position to specify the matter source which consists of the conventional 
cold matter and radiation, p = Pm + Pr with pm ~ and pr ~ a""^. Let us define 
Prad = Pr + Pdr, followiug the Conventional way that pm means the density for the cold dark 
matter as well as baryons. Then, the equation-of-state parameter uj{a) is given by 

/ \ 3r2(js ~l" f^rad 3(3, f^vac ('22") 

^ ~ 3 [^]ds + a^^^rad + + O^l^vac] ' 

where the density parameters are defined as Qi = pi/ pc evaluated at the present universe 
scale ao, which can be fixed to Oq = 1. Note that the equation-of-state parameter becomes 
a;(ao) — — ^m.n — —0.7 at the present scale, since the observational data indicates Q^s, ^rad ^ 



^^vac ^ 0.7 |26| 



In order to see how the desired evolution of the universe appears, we plot the phase 
portraits in Fig. [Uby setting fiyac — 0.7, flm — 0.3, Qk ~ 10^^/c, f2rad ~ 10~^, and Q^s ~ 



(a)k= +1 



(b)k = -l 



V 



V 



2x10-19 




2x10-19 




1x10-19 



1 X 10-19 



a 



a 



.0 



-1x10-19 



1x10-19 



-2x10-19 



2x10-19 



FIG. 1: The phase portraits for the configuration of our universe are plotted by setting Ovac — 0.7, 

nm ^ 0.3, Qk ~ 10-2 A;, Orad ~ 10-^, and fids 10-^2 for (a) k = +1 and (h) k = -1. The dots 

on the a-axis represent fixed points (a*,0) such that 1 + uj{a*) = 0. Actually, these figures looks 
the same with those of the general relativity. The crucial differences will be explicitly seen in the 
small scale region in Fig. [2l 

— 10^^^. Explicitly, for k = +1 in FiglU^a), there is a saddle node, and the corresponding 
stable and unstable manifolds are shown by thick arrows. It can be also shown that the 
universe monotonically expanding(II) and shrinking(IV) phases, and the scale factor has a 
maximum(III) if the initial expansion rate is small compared to that of the phase (II). On 
the other hand, in the region I in Fig. [TJ^a), the universe can not have a scale less than 
that of the fixed point , i.e. a > a* ^ 0.6, so that the shrinking universe starts to expand 
before the scale factor reaches a*. For A; = — 1 in Fig. [D^b), there is a center which is a kind 
of fixed points. The two straight trajectories separate closed orbits near the center from 
the curved trajectories, where the former case describes oscillating universes(II) while the 
latter one describes monotonically expanding(I) or shrinking(III) universes. Note that the 
oscillating universe is possible only when ptot < 0, and the straight line trajectories are in 
fact trivial solutions satisfying ptot = 0. As a result, the universe can experience late-time 
accelerated expansion after decelerated expansion in both cases (a) and (b), because current 
observational data indicates that the trajectory of our universe is in the region II in Fig. [T]^a) 
or I in Fig. [T]^b). 

Actually, the phase portraits in Fig. [1] are very close to the Einstein's theory in the large 
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(a)k=+l (b)k = -l 




a 



FIG. 2: The phase portraits for the configuration of our universe, r^vac — 0.7, — 0.3, ilfc ~ 
lO^^k, r^rad ~ 10~^, and J^ds ~ — 10~^^, seen near the inflation era for (a) k = +1 and (b) k = — 1. 
The dots on the a-axis represent fixed points (a*, 0) such that 1 + uj{a*) = 0. 

scale and the large expansion rate, so we have to draw small scale behaviors to expose some 
differences from the general relativity. The standard lore tells us that there are no phase 
transitions in a small scale in the Einstein's relativity, unless we consider an additional source 
such as the inflaton. However, it is possible to see a transition with the help of the dark scalar 
in the HL theory, which is shown in Fig. [2l The dark scalar looks like a quantum correction, 
in the sense that it comes from higher curvature correction and it plays a significant role 
at the small scale, a ^ 1. So, if we assume Q^s < 0, then we can have another fixed point 
near a* ~ 10~^, which gives the possibility for the first accelerated expansion. For the case 
of /c = +1 in Fig. [21 there is a center near the early accelerated era. One can now see that 
the expanding and shrinking universe in the region III in Fig. [11(a) is actually oscillating, 
which means a universe with the same energy distribution as ours should oscillate if its 
initial expansion rate is not so large enough. As for the case of = — 1 in Fig. [2l there is 
a saddle node with stable and unstable manifolds, and the two straight trajectories whose 
total energy is trivial, ptot = which separates solutions satisfying ptot > from ptot < 
solutions. 

To incorporate Fig. [H with Fig. [21 they are plotted in terms of the logarithmic scale 
in Fig. [3l and then it can be shown that the universe starts from the first accelerated 
phase followed by a decelerated expansion phase, and ends up with the second accelerated 
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(a)k= +1 

V 



(b)k = -l 

V 




ln(a) 



FIG. 3: The phase portraits for the configuration of our universe, r^vac — 0.7, Qm ~ 0.3, ~ 

IQ-'^k, ftrad ~ lO"'', and Ods 10~^^ for (a) k = +1 and (b) k = -1. The thick black (curved) 

arrows describe the trajectory of our universe, and the dots on the a-axis represent fixed points 
(a*,0) such that l+uj{a*) = 0. 

expansion. In addition to this dynamical evolution, as seen from the two fixed points in 
Fig. 131 it is also possible to obtain the neutrally stable universe described by the center at 
the early stage of universe in Fig. [3](a) and at the late time stage of universe in Fig. ^h). 
Here, 'neutral' means that a small perturbation does not decay into zero but remains small 
as time goes on. 



V. DISCUSSION 



We have studied the HL gravity coupled to the matter without the detailed balance con- 
dition in order to show the possibility to get the smooth phase transition from the first 
accelerated expansion corresponding to the early stage of the universe to the second accel- 
erated expansion throughout the intermediate decelerated expansion assuming the energy 
density for the dark scalar to be negative pds < 0. In spite of this negative energy contribu- 
tion, the total energy density ptot is positive at any cosmological scale. Note that there have 
been researches which provide concrete justifications for models with negative density, in par- 
ticular, a brane universe moving in a curved higher dimensional bulk space js?] and a model 
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of dark energy stemming from a fermionic condensate [28j. Of course, we can confirm that 
there is no such a phase transition unless the dark scalar density is assumed to be negative. 
For instance, if we take the detailed balance condition, then it simply reduces the additional 
energy contribution to p^ac = Ac'^/SvrGjv, = — Sc'^fc/STrGjva^, pdr = 27c'^fc^/327rG7vAa^, 
Pds = as seen from fundamental constants and redefinition of coefficients (jS])- In this 
case, the first accelerated expansion does not appear, that is the reason why we did not 
take the detailed balance condition. In addition, it has been shown that our universe may 
be oscillatory even with the same density parameters of the current observation, depending 
on the expansion rate. Actually, the expansion rate is related to the total energy density 
by Eq. (IT7|) . so that small expansion rate v corresponds to small h ~ ptot for a given scale 
factor. 

On the other hand, it has been well known that if there were ~ 60 e-foldings of infiationary 
expansion, then the universe was driven to be nearly fiat, removing any need for fine-tuned 
initial conditions. In addition, the infiationary expansion would have driven the density 
of magnetic monopoles to be negligible today, explaining their apparent absence. Also, 
our entire observable portion of the universe would have infiated from an initially small 
causally connected region, thereby ensuring a high degree of isotropy today. Although our 
analysis shows that the first accelerated expansion can be obtained by assuming pds < 0, 
unfortunately, it should be pointed out that it does not last for 60 e-foldings. We hope this 
problem will be discussed elsewhere. 

It is interesting to note that there exists a minimum scale amin obtained from pc > 
and Ptot > with the relation pc = ptot + Pk- Note that ptot(amin) = -pk{amm) > and 
H'^icL^in) ~ Pc (omin) = for = +1, SO the universe starts with a certain amount of energy 
and zero expansion rate, while ptot(amin) = and H'^{a^i^) ~ Pc(amin) = Pfc(a.min) > for 
k = —1, so the universe starts with zero energy density and a certain expansion rate. For any 
case, as far as we consider the positive total energy density and the critical energy density, 
there should be no initial singularity problem because of the existence of the minimum scale. 

The final comment is in order. The energy density and the pressure of the dark scalar 
which play a crucial role in this work depend on the three terms in the potential (jl]), (TiR^, 
a2RRijR^-' , asRlRjRl- One may think that three terms give rise to kinetic contributions 
so that they provide the leading spatial dependence in the gravition propagator in the UV 
region. What it means is that graviton modes may be unstable. However, this is not the case 
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since they do not modify the gravition propagator. In general, the lowest metric order in 
curvature tensors and curvature scalars is linear, which means that the three terms are not 
quadratic in metric. Of course, such terms which are quadratic in curvature, (T4Vj-RV*-R, 
a^ViRjkV'^R^'^ in Eq. (jlj), will not only add interaction but also modify the propagator; 
however, these are irrelvant to our energy density formulae f ll2p . 
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